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SUMMARY 

A Controlled Ecological Life Support System (CELSS) utilizes a plant’s natural ability to regen- 
erate air and water while being grown as a food source in a closed life support system. Current plant 
research is directed toward obtaining quantitative empirical data on the regenerative ability of each 
species of plant and the system volume and power requirements. Two techniques were adapted to 
optimize crop species selection while at the same time minimizing the system volume and power 
requirements. Each allows the level of life support supplied by the plants to be selected, as well as 
other system parameters. The first technique uses decision analysis in the form of a spreadsheet. The 
second method, which is used as a comparison with and validation of the first, utilizes standard 
design optimization techniques. Simple models of plant processes are used in the development of 
these methods. 


INTRODUCTION 

To date, life support technology is based solely on physical/chemical processes, and this is likely 
to remain true for the initial phases of the Space Exploration Initiative. However, for long-duration 
missions, such as a trip to Mars or long-term habitats on the Moon or Mars, a CELSS has the poten- 
tial to provide human life support with significant cost and safety benefits over the currently envi- 
sioned physical/chemical systems. In particular, food resupplied from Earth may be significantly 
diminished, higher plants can accomplish both air revitalization (through the release of oxygen and 
uptake of carbon dioxide) and water processing (through transpiration), and some waste disposal 
may be accomplished biologically. Figure 1 shows an example of an integrated biological and 
physical/chemical life support system (CELSS) as conceived by Dr. John Rummel, a NASA 
scientist/administrator. Other studies have resulted in variations on this conceptual design (refs. 1-3). 

Normally, a trade study is conducted to determine advantages and disadvantages of various 
design options. Trade-study techniques can be developed in parallel with research on basic perfor- 
mance parameters, so that when reliable data become available, the analysis tool is also ready to 
perform trades. These tools will become increasingly important as we begin to address the complexi- 
ties involved in integrating biological components with physical/chemical life support system 
components. 

To date, most research in the use of plants for life support has concentrated on productivity levels 
and the effects of environmental parameters on productivity. Little work has been done in evaluating 
the air and water regeneration and waste management capabilities, which would be the next logical 


81 




Figure 1. Conceptual design for a controlled ecological life support system. 


step towards developing an integrated life support system. The techniques outlined herein transform 
newly acquired plant performance data into parameters describing a CELSS system for use in trade 
studies, thus providing the link between generating data and developing an optimal CELSS design. 

We are adapting two techniques to optimize crop selection for minimum power and volume 
penalties. The first technique involves the use of decision analysis which implements a decision tree. 
The second method, intended to be both a comparison with and a validation of the first, involves the 
use of standard design optimization (linear programming) techniques. Previously, design optimiza- 
tion techniques have been applied to crop mix selection to select the minimum crop area which 
satisfies human nutritional requirements (ref. 4). While the results of that study do not account for 
power penalties, nor do they allow for air and water regeneration constraints, some comparisons can 
be made with our own data. These comparisons are reported in the results section of this paper. 

A spreadsheet is used as an interface with the user and to generate plant and system parameters. 
The user specifies the level of life support to be supplied by the plants for each life support function, 
oxygen generation, carbon dioxide uptake, water regeneration, and nutrient (carbohydrate, lipid, and 
protein) production. Plant parameters are generated both directly from input, empirical data and from 
models of plant processes. However, the techniques are intended to be generic and applicable to 
other plant parameter generation schemes. 
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APPROACH 


Derivation of Plant and System Parameters for Input to Analyses 

Parameters specifying system requirements are input into the spreadsheet section shown in 
table 1. Here the user can input which crops of those whose performance data has been entered into 
table 2 should be considered. Up to five crops can be chosen. If variety in the crop mix is desired, a 
minimum number of crops can be entered (a number greater than one will force multiple crops to be 
chosen). Currently this feature only applies to the decision analysis method. The level of life support 
to be fulfilled by plants and the daily life support requirements (ref. 5) are input into table 1, as well 
as the power and volume penalties (refs. 6 and 7) and end-to-end lighting efficiency. 

Table 2 lists parameters which were obtained from the literature and other sources of plant data 
(refs. 8-12). These parameters include each species’ rate of transpiration and biomass growth as well 
as diet composition, edible fraction, and chamber height requirement. The lighting level and photo - 
period under which these rates have been measured are also recorded and used for the power 
requirement calculations. There are a myriad of other factors that influence plant productivity, tran- 
spiration rates, edible fraction, and even diet composition, such as carbon dioxide level, oxygen 
level, nutrient solution composition, temperature, and humidity. These influences could be added in 
a more sophisticated effort to derive the parameters for the analyses input, but are not necessary for 
our purposes of technique development and demonstration. 


Table 3 lists the parameters describing plant species performance required by the optimization 
methods. Generation of these parameters can be accomplished in many different ways, from using 
empirical data to employing modeling techniques. We have elected to use a combination of these 
two methods, largely because gas exchange data for plant species are limited. Parameters which are 
more readily available from empirical data are used directly and as the basis for some simplified 
relations used to generate the remaining parameters. 


Transpiration rate and biomass production rate are taken directly from the empirical data 
recorded in table 2. The fat, carbohydrate, and protein production rates are products of the biomass 
generation rate, m bio , the edible fraction, and the fraction of the total biomass generation, which is 
fat, carbohydrate, and protein, respectively. 


fat generation 
carbo. generation 
protein generation 


= m b j 0 * (edible fraction) * 


fat fraction 
carbo. fraction 
protein fraction 


Carbon dioxide and oxygen generation rates are difficult to find in the literature, but are products 
and reactants of photosynthesis and respiration, as is biomass production. This link between biomass 
production and gas exchange rates is described by the photosynthetic equation, assuming respiration 
is ignored. The chemical reaction of photosynthesis varies with the biomass type being formed, 
whether it is fat, carbohydrate, or protein. For simplicity, here it is assumed that carbohydrate is the 
substance formed. The photosynthetic equation describing carbohydrate formation is (ref. 13). 


6CO2 + 6H2O — ^ C b Hi20 6 + 6O2 
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Table 2. Plant performance data for model input 3 
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production rales are per day 




This equation gives the ratio of moles of carbohydrate produced to the moles of carbon dioxide 
taken up and oxygen released. Converting these mole fractions to mass fractions (using the corre- 
sponding molecular weights) the oxygen generation rate, rho 2 > an d the carbon dioxide take-up rate, 
m C Q 2 , are related to the biomass production rate, m b } 0 , by 

m 02 = 1.0667 m bio 


m C o 2 =1-4 m bi0 


The power requirement is determined using the lighting level recorded in table 2 measured as 
Photosynthetic Photon Flux (PPF, the radiation given off by the lights in the wavelength band useful 
for photosynthesis), the photoperiod which is the hours each day that light is supplied to the plants, 
and the lighting system efficiency, r|, which is the end-to-end efficiency of the lighting system. The 
equation used to calculate the required power, Preq, is 

„ PPF * photoperiod 
Preq - rj 

We have not incorporated optimum lighting levels for crops in this study. For a legitimate 
application of these techniques, either data representing crops at optimal conditions must be entered 
into table 2 or a more sophisticated, compensating model must be devised. Alternatively, multiple 
entries of the same crop species could be made with parameters reflecting the crops’ performance 
when optimized for biomass production, transpiration, power, or volume conservation. 

The final calculation is the total cost. This is the sum of the power and volume (height times 
lm 2 ) times their corresponding penalties as given in table 1. 

Total cost = (PR£Q)(Power Penalty) + (Volume)( Volume Penalty) 

Note that this function could be made more sophisticated by weighting the relative importance of the 
two requirements or by adding mass as an additional cost. While mass penalties could easily be 
incorporated into the cost function, current understanding of mechanically optimized plant chamber 
mass is sufficiently limited to defer incorporating mass penalties into this study. 


DECISION ANALYSIS METHOD 

Decision analysis methods (ref. 14) provide a tool for making decisions based on a single princi- 
pal value (in our case, we have chosen to express everything in terms of cost). Most often, the tool 
used in decision analysis is a decision tree, where all possible outcomes and all possible paths to 
these outcomes are diagrammed. Many decision-tree analyses also have expected values or proba- 
bilities attached to each branch stemming from a decision or node. In our case, we are merely 
minimizing cost at each decision node, with equal probability that any particular pathway will be 
followed. The decision tree (fig. 2) is constructed such that the initial decision determines which 
crop selection is the cheapest for all permutations of solutions having the same number of crops. A 
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Figure 2. Decision tree. 


second round of decision-making is then done to determine the most economical number of crops 
one could use. 

Computation of the cost values is shown in table 4. The assumption is made that an equal area is 
allotted to each crop in a crop mix. Thus the generation/dissipation rates of the plant products, 
power, volume, and cost are averages of the individual crops in a crop mix computed on a per-m 2 
basis. The penalties incurred by this assumption are shown through the comparison of results with 
the results from the design optimization method, where this assumption is not required. The required 
generation/dissipation rates (calculated from the human requirements and the degree of support 
specified in table 1) are divided by the productivity of each crop mix to obtain the scaling factor 
(planting area in m 2 ) required to meet the specified level of life support. The largest scaling factor 
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encountered for a particular crop mix is multiplied by the cost of the crop on a per-m^ basis to obtain 
the cost penalty entered in the decision tree. 

A feature of the decision tree tool is the accessibility of cost values for all crop mixes. This allows 
the designer to investigate the cost of non-optimal solutions which might have more appeal than the 
optimal solution for qualitative reasons. For the example shown in figure 2, increasing the crop 
variety by selecting the most optimal 4 crop solution over the 3 crop solution increases the cost by 
18%. Also, if a designer preferred wheat over potatoes in the optimal crop mix of lettuce, potatoes, 
and Y (crops 1, 2, and 4), the decision tree shows an increase in cost of 14% for lettuce, wheat, and 
Y (crops 1, 3, and 4). 


DESIGN OPTIMIZATION METHOD 

An alternative to using the decision analysis method described above is to take a design 
optimization approach (refs. 15 and 16). One can then minimize the cost function while removing 
the assumption of equal crop growth areas. For example, we use a linear programming approach to 
solve the constrained optimization problem 


mini 

area 


N 


Y(co.t>, (area)i 


i=l 


Subject to the following constraints 
N 


i=l 


H 2 O transpired 
area 


N 

I 

i=l 

N 

I 

i=l 


CO 2 removed 
area 


O 2 produced 
area 


(area); > H 2 O transpiration requirement 


(area)j > CO 2 removal requirement 


(area)j > O 2 production requirement 


protein lower bound < Y f P rote ^ I l£^ ucec ! ( area ) j < protein upper bound 


i=l 


carbo. lower bound < Y f produced "j (^g^. < car b a upper bound 


i=l 
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lipid lower bound < ^ hpid ^^ UC£ --j (area); < lipid upper bound 
i=l 1 

where i identifies the crop species and N is the total number of crops being considered. The objective 
is then to identify the crop mix which minimizes cost while meeting certain requirements for air and 
water regeneration as well as food production. One could easily modify the above formulation to 
include additional requirements, such as vitamin and mineral nutritional requirements, or to include 
additional constraints, such as physical constraints on the crop-growth areas due to rack-size 
limitations or edge-effect considerations. 

This constrained optimization problem can be solved using standard linear programming 
techniques. We used a SIMPLEX algorithm, coded in FORTRAN on a MicroVax 3200 computer. 
The algorithm first identifies whether a feasible solution exists, then solves for the optimum solution 
and determines whether or not the solution is degenerate (i.e., an infinite number of solutions exist). 

The optimization method has computational advantages, especially when a large number of crop 
species or nutritional requirements are being considered. The decision analysis method requires an 
exhaustive search since the cost of each possible solution must be calculated, while the optimization 
method uses search directions to quickly find the optimum solution. Also, standard methods exist for 
performing parameter sensitivity analyses for the linear programming formulation. Such analyses 
would be very useful for performing “what if’ studies to investigate the effects of changing costs or 
productivities of the various crop species. 


RESULTS 

Table 5 shows the decision analysis output for the baseline case outlined in tables 1-3. For this 
case we have specified 100% of the requirements for oxygen, carbon dioxide, and carbohydrates to 
be fulfilled, as well as 40% of those for water, 25% of lipids, and 25% of protein. Results show three 
crops being selected (lettuce, potatoes, and Y) as the optimum mix, with a total area of 3.0 m 2 and 
total cost of 2.33. Results also show that in order to supply 100% of the carbohydrate requirement, 
oxygen, carbon dioxide, water, and protein are oversupplied, and more carbon dioxide is taken up 
than is necessary. In design optimization terminology, carbohydrates are the active constraint. 

Table 6 shows the results of the design optimization technique for the same base set of con- 
straints. Recall that this technique is not limited to equal areas for each species in the crop mix 
selected. This is reflected in the results, which in this case show three crops being selected (lettuce, 
Y, and soybeans), with a total area of 2.5 m 2 at a cost of 2.16. Examining the results, we see that the 
carbohydrate, water, and lipid requirements are active constraints. 

The optimization technique can also be used to examine the optimal solution in the case where 
protein, carbohydrate, and lipid production all become active constraints. This is accomplished by 
setting the upper and lower bounds on these variables all equal to 100% of the requirements, thus 
forcing the optimization to choose a crop selection which produces a nutritionally balanced diet (no 
over production of protein, carbohydrates, or lipids). In the baseline case, no feasible solution exists. 
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Table 5. Decision analysis baseline results 


Desfeiorr 

Analysis; 

Output 



[Total Crop 
i Number of plant . .species yfjf 
Aivaa ; plated;:; Em 23 

Total 'voiuwef/i^ 

• Total: w»t; of plants's^ 



\% Oxygen ' 
%:Carbon;: Dioxide 
% Water :V; : 

;% Lipids ' 
i % Carbohydrates 
% Protein '' 'M;:®: 


Table 6. Design optimization baseline results 


[Total cost = 2.1595 
Total area = 2.5004 m 2 



Area of crop #1 (lettuce) 
Area of crop #2 (potatoes) 
Area of crop #3 (wheat) 
Area of crop #4 (Y) 

Area of crop #5 (soybeans) 

= 0.0584 m 2 

= 0.0000 m 2 

= 0.0000 m 2 

= 2.3549 m 2 

= 0.0870 m 2 


C0 2 removed = 2.1621 

0 2 produced = 1.5724 

H 2 0 processed = 12.5040 

kg/day (216.21 

kg/day (187.19 

kg/day ( 40.00 

% of human reqt.) 
% of human reqt.) 
% of human reqt.) 

Protein produced = 

Carbohydrate produced = 

Lipid produced = 

24.48 g/day (32.21 % of human reqt.) 

633.00 g/day (100.00 % of human reqt.) 

26.25 g/day (25.00 % of human reqt.) 
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However, if crop X is substituted for crop Y a solution exists such that exactly 100% of the protein, 
carbohydrate, and lipid requirements are satisfied. In this example, lettuce, potatoes and crop X are 
selected, and carbon dioxide, oxygen, and water requirements are satisfied but are not active 
constraints. 

Note that in the base case, we have set the lipid requirement to be satisfied by the plants to a 
fairly low level (25%). In the study of reference 4, it was found that the lowest total area (only nutri- 
tional needs were evaluated) resulted when no minimum percentage of fat was made available to the 
diet via plants. We have not found the lipid requirement to be an active constraint in our case, per- 
haps because we have constrained it to only a small degree. Since additional constraints have been 
added to the baseline case (and the crops to be considered have been changed), the cost increases to 
3.0 requiring 4.8 m 2 of crop area. 

In a study by McDonnell Douglas (ref. 4), it is found that the lowest total area (only nutritional 
needs are considered) results when no lipid production requirement is placed on the plants. In the 
baseline case of our study, the lipid requirement to be satisfied by the plants is set to a low level 
(25%). For this case the lipid requirement is not an active constraint. However, when the require- 
ments for the three nutritional categories are increased to a level of 100%, lipids become the active 
constraint, which is consistent with the McDonnell Douglas study. 


SENSITIVITIES 

There are several ways in which sensitivities can be examined. Two approaches are demonstrated 
below. 

From the results of the baseline case, it was determined that carbohydrates were an active 
constraint. With this in mind, we might like to see how sensitive our answers are to the specified 
carbohydrate requirement. For example, what happens when the carbohydrate fulfillment level is 
lowered to 50%? The decision analysis results show a 40% savings in cost and a 15% decrease in 
crop area using only soybeans and Y. With this new requirements specification, the active constraint 
shifted to water fulfillment. The design optimization method showed a 38% savings in cost and a 
21% decrease in crop area. The active constraints remain carbohydrate and water fulfillment. This 
sensitivity analysis demonstrates that resupplying a portion of the carbohydrates from Earth might be 
preferred to requiring the plants to produce the entire requirement depending on resupply costs. 

Another approach to analyzing sensitivities is to examine the effect on the results when plant 
performance data for a particular species are altered. As an example, the biomass production rate and 
transpiration rate were varied (doubled and halved) about nominal values. Table 7 lists the results of 
the two methods to each of these cases . This table shows that certain increases in plant performance 
parameters result in larger cost savings than others. This method could provide a tool for steering 
plant research, as well as addressing the balance between transpiration rate and plant biomass 
production. 
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Table 7. Parameter sensitivity analysis results 


Crop 

Y 


0.35 

Biomass 

generation 

rate 

0.70 

(nominal) 

T kg 1 
Lm 2 dJ 



1.40 


Transpiration rate 
" kg ' 
_m 2 d_ 

2.5 

HBH 

10.0 

DA 

cost: 2.45 
area: 3.4 m 2 
species: LJP 

DA 

cost: 2.45 
area: 3.4 m 2 
species: L,P 

DA 

cost: 2.42 
area: 2.9 m 2 
species: L,Y 

OPT 

cost: 2.36 
area: 3.50 m 2 
species: L,S 

OPT 

cost: 2.36 
area: 3.50 m 2 
species: L,S 

OPT 

cost: 2.24 
area: 2.76 m 
species: L,Y 

DA 

cost: 2.42 
area: 3.3 m 2 
species: L,Y,S 

DA 

cost: 2.27 
area: 3.1 m 2 
species: L,Y,S 

DA 

cost: 1.87 
area: 2.3 m 2 
species: L,Y 

OPT 

cost: 2.28 
area: 2.98 m 2 
species: Y,S 

OPT 

cost: 1.96 
area: 2.34 m 2 
species: Y,S 

OPT 

cost: 1.79 
area: 2.17 m 
species: L,Y 

DA 

cost: 1.73 
area: 2.5 m 2 
species: Y,S 

DA 

cost: 1.41 
area: 2.0 m 2 
species: Y,S 

DA 

cost: 1.19 
area: 1.4 m 2 
species: Y 

OPT 

cost: 1.56 
area: 2.29 m 2 
species: Y,S 

OPT 

cost: 1.40 
area: 1.97 m 2 
species: Y,S 

OPT 

cost: 1.06 
area: 1.32 m 
species: Y,S 


2 


LEGEND: 


P 

= potato 

DA = results from decision analysis 

L 

= lettuce 

OPT = results from design optimization 

W 

= wheat 


S 

= soybeans 


Y 

= Y 
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CONCLUSIONS 


1. The design optimization method clearly provides the most optimal crop selection of the two 
methods considered. The decision analysis method results are usually within 10% of the cost and 
total crop area of the design optimization results. 

2. The advantages and disadvantages of both methods imply that the two methods are ideal com- 
panions to one another. The decision analysis method allows an interactive atmosphere with the user, 
facilitating experimentation with design specifications. A decision tree displays the cost of all crop 
mixes simultaneously. This allows the user to evaluate at a glance the cost of additional variety of 
preferred (more appetizing) crop species. 

The design optimization method has computational advantages especially when the number of 
crops being considered or the number of nutritional (or other) constraints become large. It also 
ensures that the last 10% of cost and area savings will be provided by the specified crop mix, and it 
allows the user to further constrain the optimization problem as needed. A more rigorous sensitivity 
analysis approach could also be developed using this method. 

3. There are limitations to this analysis. A significant limitation of the decision analysis method 
is the assumption of an equal area for each species in a crop mix. Both methods are limited because 
they do not allow the plant performance parameters to be simultaneously optimized (the transpiration 
rate/biomass generation rate trade-off, as well as other variables). No attempt is made to account for 
the compatibility of crop species if a common air space or nutrient delivery system is planned. Also, 
the power and volume of the processors required to support a plant system, including environment 
control, are not accounted for. 
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